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On Continuous Ambiguities in 
Mo del- Independent Partial Wave Analysis - II. 



I.N.Nikitin* 



Abstract 



A problem of density matrix determination in terms of the given an- 
gular distribution of decay products is considered. 



^ \ In the first part of this paper (Section III) a quantum system, described by 

density matrix, has been considered. In this the density matrix was diagonal 
with respect to orbital quantum numbers and the state was mixed only partially 
Now we will consider the process X — > + 0, where initial state is described by 
\ the density matrix p, for which we do not assume any other properties except 

I/"") ■ hermitity and positive definiteness. 

OS | The angular distribution is given by the expression 



CU. I(n) = J2p a pY:(n)Y p (n). (1) 

(— i \ Let us multiply both sides of this expression by Y 1 {n) and integrate this over 

the sphere 

^ 1 h = / ,Pa/3Dg0j, (2) 

r 2 

where t 7 = id n IYy =< Yy >/ are moments of distribution, 

D af3j = / d 2 n Y*YpY 7 . 



The integrals D a) 3 7 are expressed in terms of Clebsch-Gordan coefficients (see 
©§107): 



n _c_i\ra s / (2Zj + 1)(2Z 3 + 1) nh _ m2 , 

l'/|.»,.'..»>.' :1 m:, — < U \; 47r(2Z 2 + 1) *1 —mi, (37713 Zl0,i 3 
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We use spherical functions 0, different from those of [Q by a phase factor: 
= i l Y^ . Summation in (||) extends for the following values of indices: 

m 2 = mi - m 3 , l 2 = \h - h\ ■■■ h + h step 2, 

— min{Zi, I2 — ^3} < mi < min{Zi, l 2 + m^}. 
Let us represent the density matrix p a p via eigenvalues \i and eigenvectors 

i 

The form unitary matrix 

i 

The positivity of p implies that all eigenvalues are positive: A; > 0. 
Let us rewrite (||) 

h = £ M 7i X h ( 3 ) 

i 
a/3 

At 7 = from @ a normalization condition follows: 

Xi = J d 2 n 1=1, the distribution is normalized to unity. 

i 

At the given ^> l a relation (||) is a system of linear equations for A^. A matrix of 
this system is non-degenerate, except rare special values of ^ i a . 

Particularly, the matrix M 7 i degenerates when ^ = Si a . The ques- 
tion whether the matrix M 7 ; is non-degenerate reduces to clarification of 
linear independence of a set of functions ^(n)! 2 for the given complete 
set of functions \I/i(n). 

Therefore in our approach the free parameters are the unitary matrices , 
the values A.; are expressed via ^> l a from system (||), the inequalities Xi(^>) > 
select a region of possible values for <!>. 

Notion. 

The presence of continuous ambiguity in the solution of this problem is most 
obvious in n-representation of operators. The density matrix is represented by 
the kernel of integral operator 

(p^)(n) = / d 2 ri p(n,n')^(n'). 
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The distribution is equal to the diagonal elements 

I{n) = p(ft, ft) (4) 

Hence, the problem is reduced to the reconstruction of hermitean positive matrix 
from the given diagonal. Let p be a solution of this problem. The replacement 
p — > p+Ap, where Ap is arbitrary hermitean matrix with zero diagonal elements, 
conserves (Q) and hermitity of p. If the matrix elements of Ap are sufficiently 
small, then the positivity of p is also conserved. Therefore, the solution of the 
problem is not unique. 

In n-representation the analysis of p positivity and the determination of 
eigenvectors are complicated. The approach proposed above seems more conve- 
nient. 



Scheme of analysis 

When describing quantum systems via density matrix one suggests that the 
system is located in pure states with different probabilities A^. The rep- 
resentation of eigenvectors of density matrix on Argand plots is an ultimate 
goal of partial wave analysis. The input data are the moments t 1 . We propose 
to obtain the region of solution using Monte-Carlo technique. The scheme of 
analysis is the following: 

1. A random unitary matrix is generated. We use the generation of the 
columns of this matrix as a set of independent isotropically distributed 
random complex vectors with the subsequent Hilbert-Schmidt orthogo- 
nalization. 

An isotropic (Gaussian) distribution of n-dimensional vector can be ob- 
tained by generating its components as independent random values with 
Gaussian distribution: 

dp = e~ X1 ' ' dx\ ... e~ Xn2 dx n = e~ r ' 2 dV. 

The unitary matrices generated by this method are uniformly distributed 
over the unitary group invariant volume^]. 

2. The matrix M 7 j is obtained, system @ is solved for A^. 

3. The solution is accepted, if all Aj > 0, and is discarded otherwise. 

4. The eigenvalues A^ are displayed on intensity plots, the eigenvectors ^ z a - 
on Argand plots. 

The unitary transformation transfers one orthonormal basis to another. The invariant 
measure on unitary group can be defined as a product of spherical volume for one vector from 
the basis and the measure on subgroup leaving this vector invariant (little subgroup). Using 
this recurrent definition and isotropy of generated vectors, one can easily prove the above 
statement by induction. 
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Remarks. 

1) Phase rotation of eigenvectors ty l a — > e %Vi ^\ does not change the density 
matrix. One can generally choose the phases of eigenvectors in such a way that 
the component vPg in each eigenvector will be a real positive number (phase is 
measured from S-wave). Those unitary Nx N matrices are described by iV 2 — N 
independent real parameters, which being unified with N eigenvalues Aj give N 2 
degrees of freedom, contained in an arbitrary hermitean matrix p a p. 

2) Change of numeration (permutation) of eigenvalues and eigenvectors con- 
serves the density matrix. The region of solutions is symmetrical with respect 
to these permutations. Thus projections of the region onto Argand plots of 
each eigenvector coincide. One can order the eigenvalues: Ai > A 2 > ... > A^ 
and respectively permute the eigenvectors. Such ordering might lead to breaks 
of continuity in solutions. The breaks appear when a collision of eigenvalues 
occurs on a smooth solution (fig.l). The ordering transfers the continuous path 
ab to the non-continuous path ab' . One should take into account this fact when 
selecting the solutions: not only smooth paths inside ambiguity region are al- 
lowed, but also those paths, which can be smoothly continued in the points with 
coincident eigenvalues after the permutation of the correspondent eigenvectors. 

In the vicinity of the point A the eigenvectors ^i and ^2 cannot be 
close because they are orthogonal. The points A and A' cannot coincide 
on all Argand plots. Hence the break of curves on Argand plots actually 
occurs in transition through the point A along the path ab' . One can 
connect the points A and A' by continuous path, along which the density 
matrix is constant. At Ai = A2 any linear combinations of ^1 and ^2 
are also eigenvectors. The ^1 and ^2 can be permuted with the aid of 
unitary transformation, not touching others \&i (this can be developed 
continuously, because U(2) group is connective). The rotation A A' is 
developed at constant Ai = A2. Hence, though the transition aAA'b' is 
continuous, it cannot have analytical dependence on energy.) 

3) Let first N moments of distribution be measured in experiment (I < 
Lmax, N = (L max + l) 2 ). Let the unitary matrices $ have sizes N x N. We 
will solve N first equations of (0) for N unknowns Aj. The density matrices 
obtained in solution give the distributions, for which N first moments exactly 
coincide with measured ones. Other moments t 7 , 7 > N can be obtained by 
substitution of values ^ l a , Xi into (||). Generally these moments do not vanish. 
One can show from the properties of the coefficients D a py that there is a finite 
number of non-zero moments: L max < I < 2L max . 

The partial wave analysis practice is to treat higher (non-measured) moments 
as zeros with the experimental precision. One can impose conditions on these 
moments 
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where St is a statistical error of the moment. 



DY 

\SU 2 - 



N evts ' 

-^Vgvts is the total number of events, on which the moment is measured 

^ ^evts 
h = ^ U 
DKy is the dispersion of random value Y 1 (n) , 

dy 7 =< |r 7 | 2 > - | < r 7 > | 2 = y d 2 n /|y 7 | 2 - |i 7 | 2 , 
l^ 7 | 2 <-^ / d 2 ni\Y,\ 2 . 

^evts J 

Averaging both sides over the quantum number m, we have 



-r^l^ml 2 < T-Tf , i.e. |Jt, m | <-=== = Jt. 



When \tj\ > St, the deviation of i 7 from zero is statistically considerable. 

We note, that our estimate for St is I independent, all moments are measured 
with equal precision. Actually one can reliably measure a great number of 
moments (limitations are imposed only by finite angular distribution of the 
equipment: I < 2tt/ A9 Teso \). 

The contributions of the moments with high I in distribution I(n) are fast 
oscillating functions. To detect Z-th harmonics one should fill 2D histograms, 
containing at least I 2 bins. Due to limited statistics a few events will be placed in 
one bin. Hence, the presence of high harmonics can not be detected in angular 
distribution. Nevertheless, the moments of these harmonics can be precisely 
measured. 



Example. Let us consider the distribution 

I(n) = (l-a)|Yo| 2 +a|no| 2 = - a + 3a cos 2 9). (5) 



The moments are 

1 t _ a 
V47T ' 



to = / _ , t 2 a = r — , others t 7 = 0. (6) 
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Let first 9 moments t 7 be measured in an experiment (S, P and D-moments)Q 
The density matrix is described by 81 real parameters, 72 of which describe the 
unitary matrix W (\& £ 9 eigenvalues Ai can be obtained from system (||). 

Figures 2 and 3 show the result of analysis for the value a — 1/6. The 
number of random unitary matrices was 100000. In 1646 cases the solution of 
system (||) satisfied the condition A, > 0. For 54 points the additional condition 
|t im | < 10~ 2 , I > 3 held true. 

The solutions \& fill a region in the group U(9). In projection onto Argand 
plots this region maps into a set of related points. The density of points on 
figures estimates a "thickness" of the stratum, projected into the same area on 
the plot. A circle on the figures denotes the trivial solution 

Ai = 5/6 #q = 1, others #* = (S-wave) 
A 2 = 1/6 #? = 1 > others 3^ = (P-wave) 
A. ( = 0, i = 3. .9 

A few points fall in the vicinity of the trivial solution, most part of solutions is 
located in the region of close eigenvalues A, = 0...0.2. 

The points, for which the additional condition \ti m \ < 10~ 2 , I > 3 is satisfied, 
cover the same areas on the figures. A position of these solutions is analogous 
to the position of thin layer in a sphere: outer spherical layer occupies a small 
volume in a sphere, but covers the same area in projection on a plane. 



Limiting cases. 

a — -» Isotropic distribution (a — * 0) is a singular case for the given scheme 
of analysis. When i 7 = S^o/V^n, a set of equal eigenvalues A, = 1/N is a 
solution of (||) for all #. In this the density matrix is proportional to unit 
matrix: p a p = Sap/N, the isotropy of distribution (|l]) follows from the property 
l 

\ Y irn(n)\ 2 = {21 + 1)/4tt. The solution A t = 1/N has no definite limit at 

m— — l 

N — > oo and it is not of physical value. Also there are other solutions, e.g. pure 
S-wave state or 



# = 



( 1 









17(3) 
U(5) 



J 



A = 



/ Ai ... \ 

A 2 -1 3X 3 
A 3 -l 5 x 5 ■■■ 

V J 



the matrix is diagonal with respect to I (mixes m only), correspondent A; m are 
independent of to. System (|3|) has these additional solutions only if the matrix 
Mryi is degenerate. The values of parameters #, for which M 7 i degenerates, form 
in unitary group a set of zero measure. When generating matrices ^ uniformly 
distributed in U(N) volume these additional solutions will not be revealed. 

2 The distribution can differ from (jq) by the contributions of higher moments / > 3. 
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Transition to the limit a — > is shown on fig. 4. The distribution of A, 
contracts to the point A = 1/9 (N = 9), the points on Argand plots fill unit 
circles. 

The presence of such singularities indicates that the whole set of solutions 
of system (||) possesses a complicated topological structure. This structure can 
be inferred from the following examples (fig. 5): 



a) 



1 1 \ f x i \ - f l \ _ Mi \ 1 ( y 



% y J \ A2 / \ / v A2 y y — x \ ~x 



1 


1 






X 


—x 




A 2 


y 





-y J 


V a 3 




At x = y = the solutions of the system (with positive Aj) form a segment 
(0 < Ai < 1, A2 = 1 — Ai). This set has greater dimension (d = 1, a line) than 
the set of solutions at fixed x 7^ y (d — 0, a point), but less dimension than 
the set of solutions for all x,y (d = 2, a surface). Singular solutions have zero 
measure in the set of all solutions (on the hyperbolic paraboloid, fig. 5a). 



b) 



At fixed x 7^ 0, y ^ the solution is a point A^ = 1/3. When x — or y = 0, the 
dimension of the set of solutions is increased by 1, when x = y = - by 2. The 
set of singular solutions has the same dimension as the set of all non-singular 
ones (d = 2). 

In uniform random generation of parameters (a;, y) the singular solutions will 
not be found. This is essential only for singularities of the type b). 

Limits of a. Function (||) is positive on the sphere at —1/2 < a < 1. For 
these values of a the problem considered has solutions. 
Remarks. 

1) At negative a the point {Ai = 1 — a, ^1 = S-wave; A2 = a, ^2 = P-wave} 
is no longer a positively defined solution of the problem. 

2) Non-negativity of the distribution is necessary and sufficient condition 
for positively semi-defined density matrix existence, if sizes of matrices are not 
bonded (N — > 00). The necessity is obvious, the sufficience follows from the fact 
that any distribution I(n) > can be described by a pure state (see first part 
of this paper, Section I), i.e. the density matrix with single non-zero eigenvalue 



Ai = 1, Ai = 0,t > 1, *i(n) = y/lffle*™. (7) 

3) At finite N this condition is neither necessary nor sufficient. We study 
a class of distributions which can differ by harmonics with numbers greater 
than N. Even if I(n) < in some point for distribution considered, positive 
functions can exist in the same class. On the other hand, the pure states of form 
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Table 1: Number of positively defined solutions for 50000 random matrices 
with random distribution in U(9) 



a 


-1/4 


-1/8 


-1/16 





1/16 


1/8 


1/6 


1/4 




9 


2192 


14689 


50000 


14969 


2770 


823 


75 



have in general an infinite number of harmonics and are not contained in the 
considered finite dimensional class ^ l a . One can show that a set of a values, for 
which the problem has solutions in the finite dimensional class, is connective, 
i.e. it is a segment a_ < a < a+, with a+ > 1 for distribution (|J). We will not 
determine exact values a±. 

When a tends to boundary points, the volume occupied by solutions in U(N), 
tends to zero fast. The behaviour of this volume could be inferred from the data 
given in Table 1. 

At large \a\ and at a greater number of moments involved in analysis as 
well, the Monte-Carlo based approach presented here is ineffective (most part 
of solutions is discarded in positivity test Xi). One might use more advanced 
methods for solution of inequalities system Ai('S) > 0. 



Discussion 

The pure states are particular cases of mixed states, hence the ambiguity region 
for mixed states is wider than for pure ones. Even in the class of solutions, 
containing a finite number of harmonics, continuous transformations of density 
matrix are available, which do not change the distribution. One can fix the 
continuous ambiguity of solutions only imposing model restrictions on the form 
of density matrix. The illustration of this statement can be found in || . This 
work presents a formalism used in partial wave analysis of the reactions irp — > 
(3n)p, Kp — > (Kmr)p. The process is represented in a form of decay sequence^ 

meson +p—>p+ meson (J p ); 

meson (J p ) — > meson + di- meson ; 

di-meson — > 2 mesons . 

Then assumptions for the amplitude (eigenvectors of density matrix) are made, 
of which the most important are: 

3 It is assumed that the process amplitude has a general form and depends on all quantum 
numbers and kinematic variables. In this step the representation of the process as decay 
sequence is introduced for convenience and do not impose any restrictions on the form of 
amplitude. 
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1. The amplitude dependence on kinematic variables, describing the decay 
of di-meson, is introduced as a product of Brcit-Wigncr function and barrier 
factors. 

2. The amplitude of the process is represented as a product of J p state 
production amplitude and its decay amplitude. This assumption is exact only 
if a single state J p is present. 

As stated in the work, namely these assumptions fix continuous ambiguity 
of partial wave analysis (1-partially, 2-completcly). 

The author is indebted to E.B.Berdnikov, I.A.Kachaev, D.I.Ryabchikov and 
S.A.Sadovsky for helpful discussions. 

References 

[1] Landau L.D., Lifshits E.M. Theoretical Physics. V.3: Quantum Mechanics. 
Moscow, Nauka, 1974. 

[2] Aguilar-Benitez M.et al (Particle Data Group). Review of Particle Proper- 
ties// Phys.Rev. 1994. D50. p.1173. 

[3] Hansen J.D.ei al.// Nucl.Phys. 1974. B81. p.403. 

Received December 8, 1994- 



9 



Figure captions 



Fig.l. In point A the eigenvalues Ai and A2 coincide. The continuous path 
ab transfers after ordering of eigenvalues into the path ab', lying in the region 
Ai > A2 (hatched on fig.b). The eigenvectors and ^2 can not coincide (the 
case (c) is impossible). In the point A the break of curves on Argand plot occurs 
(fig.d). 

Fig. 2. Ambiguity region: eigenvalues. A circle denotes the trivial solution 
(mixture of S and P-waves). One arbitrarily selected solution is denoted by a 
star. 

Fig. 3. Ambiguity region: eigenvectors, a) b) *S? 2 - 
Fig. 4. Ambiguity region: a — > 0. 

Fig. 5. Parametric dependence of linear systems solutions (see text). 
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fig. 1 
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fig. 3a 
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fig. 3b 
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fig. 4 



